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The 56 tonne high resolution liquid scintillation calorimeter KARMEN at the beam stop neutrino 
source ISIS has been used to search for neutrino oscillations in the disappearance channel v e ^v x . 
The ^ e emitted in /x + decay at rest are detected with spectroscopic quality via the exclusive charged 
current reaction 12 C (is e ,e~ ) 12 N g . s . almost free of background. Analysis of the spectral shape of e~ 
from the i/ e -induced reaction as well as a measurement of the absolute v e flux allow to investigate 
oscillations of the type !/ e — > v T and v e — » v^. The flux independent ratio Rcc/nc of charged current 
events 12 C ( u e ,e~ ) 12 N g . s . to neutral current events 12 C ( v , v' ) 12 C* provides additional information 
in the oscillation channel v & ^v x . All three analysis methods show no evidence for oscillations. 
For the v^v T channel 90% CL limits of sin 2 (2e) < 0.338 for Am 2 > 100eV 2 /c 4 and Am 2 < 
0.77eV 2 /c 4 for maximal mixing in a simple 2 flavor oscillation formalism are derived. A complete 3 
flavor analysis of the experimental data from five years of measurement with respect to i/ e <-> v T and 
fe^v^ mixing is presented. 



I. INTRODUCTION 



One of the most interesting issues in present particle physics is to clarify the neutrino mass problem. In the minimal 
standard model of electroweak theory neutrinos are considered to be massless. But there is no compelling theoretical 
reason behind this assumption. On the other hand, in most extensions of the standard model massive neutrinos are 
allowed. A very sensitive way of probing small neutrino masses and the mixing between different neutrino flavors is 
provided by neutrino oscillations. Although there are some experimental results pointing to neutrino oscillations, it 
is still difficult to deduce a consistent and reliable set of neutrino masses and mixings 

Up to now, experimental results have been persistently interpreted in terms of neutrino oscillations using only a 2— v 
mixing scheme with a single mixing angle O. However, to perform a complete and consistent analysis of experimental 
data on neutrino oscillation search, a 3 flavor formalism of neutrinos should be adopted. In this framework, all 
experiments searching for different flavor oscillations (e.g. v^—n/e, D^-^>D e , v^-^Vr in appearance mode or i/ e —>v x , 
9 e ^D x , v^-^Vx in disappearance mode) are combined to extract global information on the oscillation parameters, 
the mixing angles as well as the mass differences Am 2 ij = \mf — m 2 \,i,j — 1, ...3. Searches for oscillations in the 
disappearance mode, although less sensitive to small mixing angles, are of great importance to restrict the allowed 
parameter space. 

In many of the 3 flavor descriptions of neutrino oscillations a so-called 'one- mass-scale dominance' Sm 2 = Am 2 12 <C 
Am 2 13, Am 2 13 « Am 2 23 = Am 2 has been adopted |3|-[|. Possible mixing to sterile neutrinos as suggested by |l(] [l^] 
is ignored whereas CP conservation is assumed, as we shall do in the following. The flavor eigenstates v a are then 
described by superpositions of the mass eigenstates Vi 

3 

v a = ^ UgjVi {a = e,n,r) (1) 

i=l 

with the mixing matrix elements U a i following the typical notation [l3| of the CKM-like mixing matrix described by 
two mixing angles The oscillation probability P a ^ for transition v a — > vp in the KARMEN experiment with 

its 'short baseline' experimental configuration (1/Am 2 k, L/E ^ 1/Sm 2 ) is then 

P aP = A- sm 2 (1.27^^) (2) 
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with Am 2 in eV 2 /c 4 , the neutrino path length L in meters and the neutrino energy E in MeV. In a simple 2 flavor 
description A = sin 2 (20) with being the mixing angle between the two v states, while in a 3 flavor analysis 
A = ^U 2 3 U 2 3 (appearance mode) or A = 4L/ 2 3 (1 — U 2 3 ) (disappearance mode) with U 2 3 = szn 2 $, U 2 3 = cos 2 ^sin 2 ^! , 
U 2 3 — cos 2 §cos 2l & ||. For an experimental oscillation search it is important to note that the second factor in (||) 
describing the spatial evolution of oscillations is unchanged in the 2 flavor analysis. In the deduction of limits on A 
and Am 2 we shall start with A — sin 2 (29) having in mind that the oscillation amplitude A has different definitions 
depending on the actual theoretical model as well as the experimental situation (appearance <-> disappearance search). 
In section |v|and[v| we shall then analyze our results in the 3 flavor scheme. 

The KArlsruhe Rutherford Medium Energy Neutrino experiment KARMEN has set stringent upper limits in the 
appearance channels jV-^e (sin 2 (26) < 8.5 • 1(T 3 for Am 2 > 100eV 2 /c 4 ) and u^v e (sin 2 (29) < 4.0 • 1CT 2 for 
Am 2 > 100eV 2 /c 4 ) [[L4j. KARMEN also searches for v e ~^v T and f e — through disappearance of v e . Here, the 
neutrinos v e emerge from /x + decay at rest (DAR) and are detected at a mean distance of 17.7 m from the source. This 
search is a complementary addition to the appearance channels. Compared to other disappearance experiments this 
search has very small, systematic uncertainties as well as an excellent signal to background ratio for the detection of 
f e -induced events. Limits on v e disappearance were obtained by reactor experiments e.g. at the nuclear power plant 
of Bugey p5[ with v e from /3-decays and the detector at 15, 40 or 95 m distance. Other limits on v e ^ v x were deduced 
from two ^"Cr source experiments in the GALLEX detector Jl(],[l7j emitting i/ e with discrete energies E v = 0.746 MeV 
(81%), 0.751 MeV (9%), 0.426 MeV (9%) and 0.431 MeV (1%). Experiments with neutrinos in the energy range of 
several GeV allow the search for appearance v^v^^ v T . In contrast to such experiments, like FNAL E531 (l8) with 
a 410 meter long decay path, in KARMEN the parameters L and E of ^-induced events can be measured with high 
accuracy on an event-by-event basis. 

The outline of the paper is as follows: We first present the experimental setup and the v e detection and identifi cation 
(section ||). Then we investigate the absolute number of identified CC reactions 12 C (v e ,e~) 12 N gs . (sect ion |III A ) 
comparing the measured cross section with theoretical predictions. In a second evaluation method (section III B[ ) we 
normalize the number of CC sequences to the number of neutral current (NC) reactions 12 C ( v , i/ ) 12 C* (15.1 MeV). 
These NC reactions are not sensitive to the v flavor. By comparing the number of CC events with NC events measured 
simultaneously systematic uncertainties in the calculation of the absolute v flux mostly cancel. A lower r atio of events 
Rexp=CC/NC than expected would also point to v e —*v x . A complementary analysis is given in section III C , where 
we search for distortions of the energy and spatial distributions of electrons from 12 C (v e ,e~ ) 12 N g s , such as would 
be due to v e — >v x . In section ^ the results of these analyses are combined and compared with other experimental 
limits for v e — > D x , v e — > v T and f„— > v e (v^^ v e ) in a 2 flavor as well as a 3 flavor mixing scheme. 



II. KARMEN EXPERIMENT 



A. Experimental Setup 

The KARMEN experiment is performed at the neutron spallation facility ISIS of the Rutherford Appleton Lab- 
oratory. The rapid cycling synchrotron provides 800 MeV protons with an average current of 200 /iA. These are 
stopped in a beam-stop Ta-D 2 0-target producing neutrons and pions. 99% of all charged pions are stopped within 
the target |l9| leading to 7r + decay at rest whereas negative pions are absorbed in the heavy target nuclei. Neutrinos 
therefore emerge from the consecutive decay sequence 7r + — > fi + + and /i + —> e + + v e + in equal intensity per 
flavor. The intrinsic contamination v e jv e from /i~ decay in the beam-stop-target is only 4.7 x 10~ 4 |19[ ). The neu- 
trino energy spectra are well defined due to the decay at rest kinematics: The from 7r + -decay are monoenergetic 
(i?, y =29.8MeV), the continuous energy distributions of v e and P M up to 52.8 MeV can be calculated using V-A theory 
(Fig. |l|a) where the characteristic Michel shape of v e from /i + DAR is given by the expression |20| ] 

2 E 
N(e)de = 4 e 2 [ 3(1 - e) + -p(4e - 3) ]de with e = — — and p = for v e (3) 

rL m ax 

The neutrino production time follows the two parabolic proton pulses of 100 ns base width and a separation of the 
maxima of 324 ns, produced with a repetition frequency of 50 Hz. The different lifetimes of pions (r = 26ns) and 
muons (r = 2.2 (is) together with this unique time structure allow a clear separation in time of the j^-burst from the 
following v e and (Fig. |l|b) which show the characteristic exponential distribution of the muon lifetime t = 2.2/xs. 
Furthermore the accelerator's extremely small duty cycle suppresses effectively any beam-uncorrelated background 
by four to five orders of magnitude depending on the time interval selected for analysis. 

The neutrinos are detected in a high resolution 56 1 liquid scintillation calorimeter segmented into 512 central 
modules with a cross section of 0.18 x 0.18m 2 and a length of 3.53m each ^T|. A massive blockhouse of 7000t of 
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steel in combination with a system of two layers of active veto counters provides shielding against beam correlated 
spallation neutron background and suppression of the hadronic component of cosmic radiation, as well as highly 
efficient identification of cosmic muons and their interactions. The central scintillation calorimeter and the inner veto 
counters are segmented by double acrylic walls with an air gap providing efficient light transport via total internal 
reflection of the scintillation light at the module walls. The event pos ition within an individual module is determined 
with a position resolution of ± 6 cm for typical energies of 30 MeV |22|] by the time difference of the PM signals at each 
end of this module. Due to the optimized optical properties of the organic liquid scintillator p3| and an active volume 
of 96% for the calorimeter, an energy resolution of <je = 11-5%/ y/ E[M eV] is achieved. The KARMEN electronics is 
synchronized to the ISIS proton pulses to an accuracy of better than ±2 ns, so that the time structure of the neutrinos 
can be exploited in full detail. 

B. Detection and Identification of v c 

Electron neutrinos v e from fi + decay at rest (DAR) are detected via the exclusive charged current (CC) reaction 
12 C ( v e , e~ ) 12 N g . s . which leads to a delayed coincidence signature. The ground state of 12 N decays at its production 
point with a lifetime of r — 15.9 ms: 12 N g s — > 12 C + e + + v e . The detection signature consists of a prompt electron 
within a few /is after beam-on-target (0.6 < t pr < 9.6 (is) followed by a spatially correlated positron (e~/e + sequence 
with 0.5 < tdiff < 36ms). To reduce cosmic background, no event is accepted within a software deadtime of 20 /is 
after any signal in the detector system. Remaining sequences induced by cosmic background are studied with high 
statistical accuracy by demanding the prompt event to emerge in a 200 /is long time interval before beam-on-target 
when no neutrinos are produced at ISIS. This CC reaction has been previously studied in detail |p4| , with respect to 
weak nuclear form factors |2{J and to the CC helicity structure of muon decay | j26| . 

The data set used for this v e ^v x analysis is taken from 1990-1995 corresponding to 9122 C protons on the ISIS 
beam stop target, or 2.51 x 10 21 /i + decays at rest in the target. Figure || shows the signatures of 513 events surviving 
all cuts. These events contain 499.7 ± 22.7 ^-induced e~/e + sequences and 13.3 ± 0.8 cosmic induced background 
sequences. The energy spectra of the prompt e~ and the delayed e + follow the Monte Carlo (MC) expectations 
from fi + DAR and the 12 N /3-decay. The time distribution of the prompt event reflects the muon lifetime, the time 
difference tdiff shows the 12 N lifetime folded with the detection efficiency. 

The stringent signatures of the 12 C ( v e , e~ ) 12 N g s . reaction lead to an excellent signal-to-background ratio of 37:1 
which allows identification of v e almost free of background. Accordingly, the search for v e disappearance has very 
small systematic errors and the sensitivity in this oscillation channel is not limited by background events. 

III. OSCILLATION SEARCHES 

Flavor oscillations of the type v e -^ or v e -^ v T will result in a smaller number of detected 12 C ( v e , e~ ) 12 N g s . se- 
quences than expected because the analogous charged current reactions 12 C ( , ^~ ) 12 N g . s . and 12 C ( v T , r~ ) 12 N g . s . 
cannot be induced due to the small neutrino energies (E u < 52.8 MeV). In addition to the comparison of the observed 
v e rate with expectation, the effects of v e disappearance can be investigated by a detailed study of the remaining 
^-induced events. First, the characteristic energy dependence of flavor oscillations (Q) will lead to a distortion of 
the electron energy distribution. Secondly, the 1/r 2 dependence of the v e flux from /i + DAR in the beam stop target 
will be distorted. The specific features of neutrino production at ISIS and detection by KARMEN require the two 
oscillation channels v e ^v^ and v e ^v T to be discussed separately. 

We first consider transitions v e ~^v T . As described above, v T cannot be detected via CC reactions, so oscillations 
v e — >v T directly lead to a reduction of CC sequences and consequently of the measured flux averaged cross section 

(fCC )cxp- 

The second channel, v e ^ mixing, is more complex due to the simultaneous presence of in the neutrino 
beam from /i + decay. Assuming CPT invariance, there would be a second source for CC reactions in the KARMEN 
detector, namely v e from D^— > D e leading to 12 C ( D e , e + ) 12 B since from /i + DAR at ISIS are produced with the 
same intensity and at the same time as v e . These would mostly compensate the reduction through v e — tv^ of the 
measured cross section due to the following reasons: 

Because 12 B belongs to the isobar triplet 12 B- 12 C- 12 N, the Q-value for 12 C ( v e , e+ ) 12 B (Q = -13.4MeV) as 
well as the 12 B lifetime (r = 29.14ms) and the endpoint energy of 12 B g . s . — > 12 C + e~ + D e (E — 13.4MeV) are 
very similar to the corresponding quantities for 12 C ( v e , e~ ) 12 N g . s . with the subsequent 12 N /3-decay (see Fig. ||). 
Moreover, the cross sections are similar with (<r) = 9.0 • 10~ 42 cm 2 [p7| induced by D e from P^—tPe with large Am 2 
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(see Fig. ^). So, as one CC reaction disappears, the other takes its place. [] Figure |^ demonstrates this argument 
in detail showing the relative detection efficiencies e(y) and e(P) as a function of the oscillation parameter Am 2 . 
Taking into account the different energy dependent cross sections for both CC reactions, the detection efficiency 
e(u,Am 2 ) within the KARMEN detector for 

Ve~ * v x disappearance into v T or is only slightly higher than for 
v^—tVe appearance e(u,Am 2 ) which also incorporates the different detection efficiency for the sequential e + /e~ due 
to the different lifetimes and Q-values of 12 N/ 12 B. The different v energy spectra and energy dependences of the 
cross sections are responsible for a small shift in the position of minima and maxima of the efficiencies for v e — > v x and 
v^— > v e . Therefore, at Am 2 m 5eV 2 /c 4 a small net increase of CC reactions is expected assuming mixing between v e 
and v^. 

Following these considerations the performed search for v e — > v x must actually be interpreted in terms of z^ e <-> and 
v e <-» v T mixi ng where the z^ e — >z/ r search is much more sensitive ( in terms of the absolute number of CC events; see 
section III A) than v e -^ due to the almost complete cancellation by the mirror v^— > i> e . An analysis of the spectral 
shape of CC sequences, however, is sensitive to both mixings, v e ^v T and z/ e <-»i^, as will be shown in section III C . 



This shape analysis is performed in a complete 3 flavor scheme where the amplitudes for the mixings v e ^> v T and 
Ve^Vy, are completely free and independent parameters. 

A. Absolute Cross Section 

The comparison of the number of ^-induced CC reactions with the expectation from theoretical calculations is 
done on the basis of the measured flux averaged cross section ( ctcc )oxp for v e from /i + DAR. Taking into account 
the overall detection efficiency of tec = 0.328 the 499.7±22.7 z^ e -induced events (see Fig. ||) correspond to a flux 
averaged cross section of 

( <tcc )cx P = (9.4 ± 0.4 (stat .) ± 0.8 (ss/st .)) x l(T 42 cm 2 (4) 

where the systematic error is almost entirely due to the uncertainty in the calculation of the absolute v flux from 
ISIS |l9| . The measured cross section is in excellent agreement with theoretical calculations based on different models 
to describe the 12 C nucleus (see table |), i.e. there is no indication of v e disappearance. For quantitative calculations 



we use the mean theoretical value of (o'cc)th from |27-31| with a realistic estimate of the systematic error pS 



[ ccc )th = (9.2 ± 0.5) x 10 42 cm 2 . The v e disappearance oscillation probability Pv e — > v x can be written in terms of 
cross sections as 

p _ , { O-CC )exp (t ., 

n/ e -> v x — J- - -j— r — \o) 

\&CC )th 

1 9 - 4±0 - 89 M 

= 1 -^2±0J (6) 

where the statistical and systematical errors of the experimental value have been added quadratically. To get a density 
distribution for values of Pi/ e —n/ x we sampled (acc )cxp and (ace }th from Gaussian distributions with the given 
mean values and widths (Fig. ^). The resulting distribution of P Vfr + Vx is slightly non Gaussian due to the ratio 
introduced in (^|). Only positive values of Pv e —n/ x correspond to the physically allowed region of oscillation. To 
extract a 90% CL upper limit we therefore renormalize the distribution in Py e -^> v x following the most conservative 
Bayesian approach The limit deduced is then 

Pv e ^v x < 0.169 (90% CL) (7) 

The interpretation of this limit in terms of the oscillation parameters Am 2 and mixing angle is given in section ^vj . 

B. Ratio of Cross Sections CC/NC 

KARMEN measures, simultaneously, the CC reaction 12 C (is e , e~ ) 12 N g s . with its sequential e~/e + signature and 
the NC reaction 12 C( v , v 1 ) 12 C* ( 1 + , 1 ; 15.1 MeV ) induced by v e , from /i + DAR. The NC event signature consists 



1 The detection of — > v e oscillations by p ( v e , e + ) n leads to a different signature as the neutron capture time t„ = 120 /is 
is much faster than the /3-decays of 12 N or 12 B (0.5-36 ms time interval). 
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of a 15.1 MeV 7 from deexcitation of 12 C* detected within the v e: v u time window after beam-on-target. This reaction 
was observed by KARMEN for the first time |3(| and is described in detail in |57],[38| . The updated cross section for 
the 1990-1995 data 



( ctnc }cx P = (10.8 ± 0.9 (rtot0 ± O.S (sys t.)) x 10" 



■ 42 cm 2 



(8) 



again is in very good agreement with theoretical calculations (table |l|) . 

In the case of oscillations between v e and other v flavors the NC cross section remains unchanged due to the flavor 
universality of NC interactions. The number of CC events 12 C ( v e , e~ ) 12 N g s . would decrease due to v e ^ v T . By 
looking at the ratio of observed events (i.e. the measured cross sections) the systematic uncertainties in ( a } CX p are 
significantly reduced so that statistical fluctuations dominate the error on R exp : 



Rexp — 



cxp 



(CNC ) 



: 0.86 ±0.08 



(9) 



cxp 



This v flux independent ratio R exp is in good agreement with theoretical predictions (see table |). Again, there is no 
indication of v e disappearance. 

As theoretical prediction for R t h we combine the -R-values from [p7|-j3l"l in a straight forward way and take the 
mean value and its variance as an estimate for the systematic error: R t h = 0.91 ± 0.035. The v e ~^ v x oscillation 
probability can then be written as 



1 — Rexp/ Rth 

I 0.86 ±0.08 
~ 0.91 ±0.035 



(10) 
(11) 



In close analogy to the evaluation in section III A an upper limit can be extracted after renormalization to the physical 



region: 



< 0.192 



(90% CL) 



(12) 



This limit is slightly higher than the limit Mj in section III A due to the relatively high value of the experimental 
cross section (ctnc )exp compared to the theoretical predictions. 



C. Spectral Shapes 

The third means of searching for oscillations v e — > v x is to carry out a maximum likelihood (ML) analysis of the 
measured energy and spatial distribution of the prompt electrons from 12 C (u e) e~ ) 12 N g s . reactions. This method 
is not sensitive to the absolute number of detected sequences but to distortions of the expected spectra in energy 
and position within the detector volume. In order to minimize contributions from cosmic induced background we 
applied more stringent cuts in time (corresponding to 3 ) and spatial distribution (restricted fiducial volume of 
86% of the central detector) for the prompt event resulting in 458 accepted sequences (see Fig. [7]) with 5.8 background 
events determined in the appropriate prebeam evaluation. With a neutrino signal-to-background ratio of about 80:1 
the cosmic background can be included as a fixed component in the likelihood analysis and will not be discussed 
further. The aim of the ML analysis is to calculate a combined likelihood for all 458 events on the basis that the 
two-dimensional density distribution f(L,E ' ) is sensitive to varying amounts of oscillation events from v e —*v Tl 
Ve-^v^ and D^^>D e . E & — is related to neutrino energy by E - = E Ve — 17.3 MeV, = Ly e since the electron is 
detected at the interaction point. The absolute energy scale is known to a precision of ±0.25 MeV from analysis of 
Michel electrons from decay of stopped cosmic muons. The density function / for no oscillations is 

f noose. (L, E e - ) = -J— • n(L) ■ ® Ve {E Ve ) ■ a(E Ve ) (13) 

"noose. 

with the normalization factor F integrated over the v e energy and the detector volume 



2 Note that this argument is only valid for flavor oscillations with neutrinos v ei v^,v T . Neutrino oscillations v e ~^ v s involving 
sterile neutrinos would reduce both the number of i^-induced CC as well as v e ,Vp, -induced NC reactions. 
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F noosc .= [ [ n(L) ■ ®v e {E Ve ) ■ <7{E Ve )dLdE Ve (14) 

JEy JDet 

The function n(L) represents the spatial distribution of events within the rectangular detection volume resulting from 
an isotropic neutrino flux from ISIS ~ L~ 2 ). Q Ve denotes the v e energy distribution from /i + DAR (see Fig. |l]a), 
a{Ep e ) the energy dependent cross section for 12 C (u e ,e~ ) 12 N g . s . (see Fig. [|). 

For an oscillation v e — > v x with given Am 2 the density function f Am 2 normalized by the integral F% 2 over L and 
E can be written as 

r An AL,E e -) = -i— ■ «n 2 (1.27^^) • n(L) ■ *„ (25„ e ) ■ a{Ey e ) (15) 
The appearance of v e through v^— > v e is described by the density 

flm^L, E e+ ) = -J- • Sln 2 (l.27^^) • n(L) • SptepJ • aLEpJ (16) 



in close analogy to (|15[). The density functions used in the ML analysis are MC-simulated electron/positron (see 
Fig. ||) an d positron/electron (see Fig. ||) distributions including all effects due to detector resolutions and the cuts 
applied . The combined likelihood for the N = 458 events is calculated with respect to the fraction r*i of initially 
expected electrons now absent due to v e — > v T and r 2 of the combined transition v e — > and v^— > ^ e taking into account 
the different detection efficiencies e{y, Am 2 ) and e(p, Am 2 ) (see Fig. ||): 



AT JV 



Wtn,^) = ]J L k Am2 {r 1 ,r 2 ) = JJ { (1 + ri + r 2 ) • fnoosc.(E k , L k ) 

k=l k=: 

~ r i ' fl m ^ E k,L k ) - r 2 



k=l k=l 

1 



(z/, Am 2 ) — e(^, Am 2 ) 
■ [e(^Am 2 )-/X m2 (i? fc ,L fc )- e (P,Am 2 )./X m2 (i? fe ,L fc )]} (17) 

In this specific construction of the likelihood function L the event number is fixed at N = 458, only the shape 
differences of the components are analyzed. Although not explicitly quoted in ( |l7j) we also include into L^m 2 ( r ii r 2) 
the knowledge, that the delayed energy E^ci and the time difference idiff differ for 12 C ( v e , e + ) 12 B (see Fig. |^). 

Maximising £Am 2 ( r ii r 2) gives the most likely oscillation fractions ri max (Am 2 ), r2 maa; (Am 2 ). This optimization is 
done separately for each value of Am 2 in the range 0.1 to 100eV 2 /c 4 . Over this range of Am 2 the fully 2-dimensional 
best fit values are compatible with n = r 2 = within the ltr error band (see Fig. |^). As there is no indication 
for oscillations, upper limits n(90%CL), r 2 (90% CL) for each Am 2 are determined. In the case of v e ^v T the ML 
analysis is not able to discriminate oscillation events for 2 < Am 2 < 3eV 2 /c 4 and Am 2 > 50eV 2 /c 4 . This is due to 
the fact that energy and spatial distributions for these oscillation cases are almost identical with the expected electron 
spectra without oscillations. However, for v e <^>Vn mixing, the shape analysis is also sensitive for large Am 2 since we 
expect events with E pr > 36MeV from P^^Pe-IJ The 90% CL limits for each Am 2 are calculated by integrating the 
corresponding 2-dimensional likelihood function for r\(r 2 ) over the whole r 2 {r\) range and then renormalizing to the 
physical region r\(r 2 ) > applying a Bayesian approach near the physical boundary r\(r 2 ) = 0, respectively. The 
limits obtained are equivalent to the limits Pv e —n/ T (90% CL), P Ve -^ v (90% CL) in oscillation probability. The 
translation of these limits into 90% CL limits for the neutrino oscillation amplitude and mixing angles for different 
is shown in section 

IV. EXCLUSION LIMITS 



In this section, we translate the results of the different methods obtained in section III as limits on the oscillation 
probability (cross section evaluation) or on the j/-mixings (shape analysis) into limits of the oscillation amplitude. 
This oscillation amplitude A is then expressed in terms of the mixing angles, for example A eT = sin 2 (20) (2 flavors) 
or A eT — Asin 2 <&cos 2 <&cos 2 ^ (3 flavors, see section |). 



3 The energy window 36 < E pr < 50MeV can in fact be used to search for — > v & oscillations in the appearance mode with 
? e detection via 12 C (v e ,e + ) 12 B. 
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A. Exclusion limits from measurement of the CC cross section 



We first discuss the limit obtained from the evaluation of the absolute number of detected CC sequences (section 
ilIA|) : iy.-n,, < 0.169 (90% CX)0 in a 2 flavor as well as in a complete 3 flavor scheme. 



For a fixed neutrino energy and a fixed source-target distance L, the oscillation probability P is related to the 
mixing amplitude A via the spatial evolution term in (||). For an experimental configuration with continuous v 
energies and a large volume detector with its resolution functions, this spatial evolution term is determined by MC 
simulations resulting in oscillation detection efficiencies e(v, Am 2 ) and e(z7, Am 2 ) as shown in Fig. |5|. The upper limit 
of the oscillation amplitude for a given Am 2 is then 

A eT ■ e(y, Am 2 ) < P Ve ^ Vt (90% CL) (18) 

For Vp-^v^ the analogous MC simulations of positrons from P^—tPe based on /^ m2 were performed resulting in 
e(P, Am 2 ) so that 

A^ ■ (e(u, Am 2 ) - e(P, Am 2 ) ) < Pj, e -» V/t (90% CL) (19) 
Pv e —>i/ T (90% CL) and Pp e ^ v (90% CL) are the evaluated limits of the oscillation probability from the different 



analyses of section [II 



Trivially, Pv e — > u x — Pv e -^> v T +Pv e -^ Vn assuming three neutrino flavors. However, in a 2 flavor scheme ignoring a 
second mixture, Pv e —> v x = Pv e —> v T or Pv e — * v x = Pv e -^> Un respectively. Therefore the 90% CL limits in sin 2 (29) 
can be evaluated following @ and ©: sin 2 (29) < 0.169/e(i/, Am 2 ) for v e -+v r or sin 2 (2G) < 0.169/(e(i/, Am 2 ) - 
e(P, Am 2 )) for v e —^i/^. Figure [To| a shows the exclusion curves from this evaluation. Note that logical consistency 
requires that either one or the other limit holds, not both simultaneously. Parameter combinations of Am 2 and 
sin 2 (20) to the right of the curves are excluded at 90% CL. Taking v e ^ v T mixing, this corresponds to mixing angles 
of sin 2 (20) > 0.338 excluded for large Am 2 and Am 2 > 0.77eV 2 /c 4 excluded for maximal mixing. Due to the 
very small efficiency e(v, Am 2 ) — e(p, Am 2 ) only three minor parameter regions can be excluded assuming v e -^v^ 
mixing. Figure |To| b shows the KARMEN exclusion (curve 1) for v e ^ v T together with exclusion limits from other 
experiments. Curve 2 represents the 90% exclusion limit from a negative oscillation search in the channel P e — ► x at 
the nuclear power reactor of Bugey, France pl| . The second limit of a disappearance search comes from the G ALLEX 
detector in which the v e flux from the two 51 Cr source tests has been investigatecj^] with respect to v e -^v x . 
The Fermilab exp eriment E531 searched for CC r-production through is e —>- and v e — * v T appearance in a wideband 
v^(v e ) beam |18|| . The exclusion limit from the v e ^v T search is shown as curve 4. Due to the small experimental 
values L/E < 0.1 with a most probable value of L/E « 0.01, E531 is only sensitive to values of Am 2 > 10eV 2 /c 4 . 

In a more detailed 3 flavor evaluation, the KARMEN oscillation limit Pv e -^> v x corresponds to 

A eT e(v, Am 2 ) + A^ (e(i/, Am 2 ) - e(u, Am 2 )) = P Ve ^ Vx < P Ve ^ Vx (90% CL) = 0.169 (20) 

In Fig. [ll] the correlation of A eT and A efJ - limits is shown for some examples of Am 2 . The parameter A e ^ is 
constrained only for regions at Am 2 w 2-3eV 2 /c 4 , Am 2 w 7-9 eV 2 /c 4 and at Am 2 w 13 eV 2 /c 4 (see Fig. |Io|a), whereas 
for Am 2 «5eV 2 /c 4 the upper limit of A eT even increases with A efl due to the fact that e(c,Am 2 ) — e(i>, Am 2 ) is 
negative. 

From A eT = Asin 2 <frcos 2 <frcos 2 ty we can evaluate the 90% CL limits for the mixing angles $ and \& taking the 
largest upper bound of A eT with respect to A e ^ . With the double logarithmic presentation of tan 2 tan 2 ^ we 
follow again H. Figure |l2|a,c shows the 90% CL limits in A eT and A e ^ from the cross section analysis as dashed 
curves (labeled with a). In Fig. |l2]b, these limits from A eT are plotted as dashed curves for two examples of Am 2 
(2eV 2 /c 4 and 100eV 2 /c 4 ) excluding areas in ($, \f) to the left. The corresponding limit on ((I",^) from A efJ/ for 
Am 2 = 2eV 2 /c 4 is shown in Fig. [T^ d (dashed curve labeled a). 

The upper limit for v e ^v T: P Ve -^ v < 0.169 (90% CL), is reliable due to the spectroscopic measurement of the 
i/ e -flux from u + DAR at ISIS and the nearly background-free 12 C (v e ,e~ ) 12 N g s . detection reaction. This limit is in 
contradiction to a model of 3 flavor neutrino mixing and masses calculated by Conforto |39|,^0|. Although being in 



4 Here, 
above. 



we do not consider the limit (113) from section 



IIIB 



since it is less stringent and therefore included in the limit used 



The complete analysis of the two Cr source tests give s a v signal of R = 0.93 ± 0.08 compared with the expectation and 
results in slightly higher upper limits than an analysis |l7| of a preliminary result of the first source test with R = 1.04 ±0.12. 
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conflict on a three standard deviations level with the negative result in v e ^ v T of FNAL E531 §J this compilation 
of different experimental results predicts a transition probability Pv e -^> v T = 0.23 ± 0.06 (90% CL) expected in the 
appropriate 'short-baseline' regime. Figure |l|a shows the calculated point {Am 2 ,A eT ) = (377±29eV 2 /c 4 ,0.46±0.12) 
lying completely in the area excluded at 90% CL by the KARMEN analysis of the absolute number of CC events. 



B. Exclusion limits from the spectral shape analysis 



In the 2-dimensional complete shape analysis (section III C ) the oscillation probabilities P Ve -^, v and Py e ^y T 
themselves depend on Am 2 . The limit in the mixing amplitude is then 

A eT < ^%CL,Ar^) ^ ^ M90% CL A m p 

e(v,Am 2 ) e(is,Am 2 )-e(D,Am 2 ) y ' 



The resulting limits in A eT and A e ^ are shown in Figure |12p , and Figure |12p , respectively, as solid exclusion curves 
(labeled ML) together with the exclusion curves evaluated above (dashed lines). For v e -^ v r the analysis of the spectral 
shape of L and E is only sensitive and competitive to the evaluation of the absolute event number in regions of about 
0.5 E/L < Am 2 < 5E/L corresponding to Am 2 values of « 3 to 30eV 2 /c 4 for this experiment whereas for v^v^ the 
shape analysis is far more sensitive. This is due to expected events from v^— > D e with positron energies E e + > 36MeV 
(see also Fig. ||), where electrons from 12 C ( v e , e~ ) 12 N g s cannot arise. Fig. |l2| b and Fig. p~2] d show the regions 
in (<&, "J 7 ) space for various Am 2 values and Am 2 =2eV 2 /c 4 , respectively, excluded by measurement of v e ^v T and 
Ve^+ mixing by KARMEN. These limits are compared with the exclusion limits from the Bugey v e disappearance 
experiment where the particular value of Am 2 chosen for comparison, 2eV 2 /c 4 , has no special significance. Note that 
there are no limits on A eT from accelerator experiments like FNAL E531 below Am 2 = 10eV 2 /c 4 . 



C. Global limits from the 3 flavor analysis 



Up to now, there are no bounds from v e ^v T to the mixing angles <f>, \& in the region Am 2 < 10eV 2 /c 4 (see for 
example FNAL E531 limit in Fig. ^djp). With A eT = 4U 2 3 U 2 3 = 4sin 2 $cos 2 §cos 2l S> (see equ. ||) we can set new 
limits to the mixing angles <!> and & from v e -^ v T oscillation search (see Fig. [12|d, [l3|). For a given Am 2 of for example 
Am 2 = 2eV 2 /c 4 the limits on the oscillation amplitudes A eT and A e ^ from the v e — >v x search can be combined 
with other limits from KARMEN as well as from other experiments. Figure |l3| demonstrates the complementarity of 
different experimental results with respect to the mixing angles. 

The shaded area represents the 90% likelihood regions obtained by LSND [|d| in favour of v^— ► v e oscillations. 
The other curves show 90% CL limits from KARMEN (D u — » V e appearance search, v e — and v e — >u T from this v e 
disappearance search), BNL E776 iy^— > v e appearance f42[| ) and FNAL E531 (^ M — > v T appearance |l8| ). For this value 
of Am 2 , the — > v e exclusion curve from BNL E776 just covers the allowed region by LSND[] with the KARMEN 
v^^D e sensitivity now substantially improved following an experimental upgrade (see section ^). However, the 
complementary searches of KARMEN v e -^r v T and FNAL E531 — > v T clearly exclude regions of the LSND evidence 
with tan 2 $> ~ 1 and tan 2 ^> « 1, respectively. Again, the necessity of investigating all oscillation channels to aim for 
a consistent description of neutrino mixing becomes evident. 



V. CONCLUSION AND OUTLOOK 



The KARMEN experiment has found no positive evidence for v oscillations in the investigated channels v e ^> v T 
and v e —*v^ through disappearance of v e . Charged current reactions induced by v e were analyzed by three different 
methods investigating the absolute number of events, the ratio of CC to NC events and the spectral shape of energy E 
and target distance L. The sensitivity of the v e — > v x search is limited by the accuracy of the theoretical cross section of 
12 C ( v e , e~ ) 12 N g . s . and the v flux calculation in the first evaluation, and by the relatively small statistics of 458 events 
in the case of the shape analysis. The experimental parameter L/ E can be measured with high accuracy on an event 



6 Note that for values Am 2 < 2eV 2 /c 4 BNL E776 with its 90% CL limit does not cover the whole LSND region in ($,*) any 
more. 
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basis and differs significantly compared to other v e disappearance searches at reactors or v e — > v T appearance searches 
at accelerators. In the maximum likelihood analysis, the 3 flavor approach assures fully independent treatment of 
i/ e <-> v T and v e ^ mixing as free parameters. The sensitivity of KARMEN is competitive with existing limits from 
oscillation searches at reactors and improves limits from accelerator experiments for v e ^ v T for Am 2 < 10eV 2 /c 4 . 

In 1996, the KARMEN detector system was improved by an additional veto counter system consisting of 300 m 2 of 
plastic scintillator modules placed within the massive iron shielding of the central detector. This new veto system will 
reduce significantly cosmic induced background for the oscillation search, especially in the appearance channel Vu,— > v e 
fhjf . In the ongoing measuring period we expect again about 400 i/ e -induced CC sequences from 12 C (v e , e~~ ) 12 N g s .. 
Therefore, the statistical significance of this v e —> v x investigation will slightly improve ending up with about 900 
i/ e -induced CC reactions in the KARMEN detector after a measuring period of 2-3 years. 

We gratefully acknowledge the financial support of the German Bundesministerium fur Bildung, Wissenschaft, 
Forschung und Technologie (BMBF), the Particle and Astronomy Research Council (PPARC) and the Central Lab- 
oratory of the Research Council (CLRC) in the UK. 
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model 


12 C(*/ e ,e-) 12 N g . s . 




R 


reference 


EPT' 


9.19 


9.87 


0.93 


Fukugita 




EPT 


9.0 


10.6 


0.85 


Mintz (28 




CRPA S 


9.29 


10.53 


0.88 


Kolbe §9 




Shell with OBD a 


9.1 


9.8 


0.93 


Engel [BC 




Shell with OBD lu 


9.4 


9.9 


0.95 


O'Connell 131) 


EPT 




10.3 




Bernabcu J'i.'.f 


Shell 




9.9 




Parthasarathy [34|] 




9.4 ±0.4 ±0.8 


10.8 ±0.9 ±0.8 


0.86 ±0.08 


this experiment 



TABLE I. Comparison of theoretical cross sections (er)th in units of 10 _42 cm 2 averaged over the neutrino energy from 
DAR and the values (ct) cxp measured with KARMEN (including statistical and systematical errors); for the NC cross section 
u—Ue+u^ in case of no oscillations; R indicates the ratio acc ( v £ )/onc (v e + v^i ) ■ 
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FIG. 1. Neutrino energy spectra (a) and production time (b) at ISIS. The proton double pulses are repeated with a frequency 
of 50 Hz. 



7 Elementary Particle Treatment 

8 Continuum Random Phase Approximation 

9 One Body Density 

10 This value for the CC exclusive cross section is calculated by |3(]] using the OBD of |HJ], but is also obtained with the 
original computer code NUEE from |32j. 
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FIG. 2. Signatures of the v £ detection reaction 12 C ( v e , e~ ) 12 N g . a . : time and energy distributions of the prompt electron 
(a;c) and the delayed positron (b;d). Data points are with background subtracted, MC expectations are shown as histograms, 
total background as shaded histograms. 
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FIG. 3. Simulated detector response in energy, target distance and time difference including readout efficiency of 
12 C ( v £ , e~ ) 12 Ng. s . (solid histograms) and 12 C ( v e , e + ) 12 B (dashed histograms) for Am 2 = 100 eV 2 /c 4 normalized to the same 
arbitrary number of events. 
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FIG. 4. Energy dependence of the cross sections of charged current (CC) reactions 1 C(u e ,e ) 1 N g . s . and C ( v e , e + ) 1 B 
and of neutral current (NC) 12 C ( v , v' ) 12 C* ( 1+ , 1 ; 15.1 MeV ) for different v flavors in an elementary particle treatment of 
the Carbon nucleus P71| . 
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FIG. 5. Relative detection efficiencies e as functions of the mass parameter Am 2 : a) e(v) for v e ^v T mixing; b) v^^v^ 
mixing: curve 1 shows t(y) for v e —*Vn disappearance, curve 2 shows e(u) for v^^Ue appearance; curve 1-2 shows the net 
efficiency e(v) — t{v) for the combination of and v^—* v e - The efficiencies e are normalized to the detection efficiency of 

12 C ( v e , e~ ) 12 Ng. s . sequences within the applied cuts. 
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FIG. 6. Simulated distributions of cross sections (o"cc)cx P and (crcc)th of 12 C ( v e , e~ ) 12 N g . s . and resulting Pv e ^u x 
according to Gaussian distributed values of (ucc }. The shaded areas show the physical region, darkly shaded are the upmost 
10% of the shaded P v ^ Vx distribution. 
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FIG. 7. Distributions of energy and target distance of prompt events (data points) with stringent cuts for the v e — ► u x shape 
analysis and underlying MC expectations without oscillations. 
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FIG. 8. Examples of simulated e~ energy and distance distributions (E— and L-projection of f noose. —/Xm 2 ) m ^ ne oscillation 
scenario v e ^>v T for maximal mixing A eT = 1 with different Am 2 as well as for no oscillation. The spectra show the visible 
electron signal in the detector and are normalized to the experimental event number iV = 458 for no oscillation (see Fig. . 
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FIG. 9. Result of the 2-dimensional ML shape analysis of the 458 events in the channels v T and f M + P e . 
Shown are the best fit values for oscillation contributions (solid lines) with the corresponding 1-a errors (dashed lines) as a 
function of the oscillation parameter Am 2 . 
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FIG. 10. 2 flavor sensitivity plot for the different channels Ve^v^ and f e — > v T from KARMEN (a) and comparison with 
limits of other experiments (b). 

a: KARMEN limits from the absolute cross section for v R ^> v T (solid line) and v^-^Vy, (dashed line) mixing, 
b: Limits from KARMEN u e -n> T (l), Bugey u e -*D x (2), GALLEX u e -n/ x (3) and FNAL E531 v e -^u T (4). 
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FIG. 11. 3 flavor 90% CL exclusion limits of the mixing amplitudes y4 eM and A eT based on the extracted oscillation limit of 
Pv e ^ fx < 0.169 (90% CL) from the cross section analysis (see equ. M). Areas to the right of the lines are excluded. 
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FIG. 12. a,c: KARMEN exclusion limits of the mixing amplitude A eT and A eM (solid line = shape ML; dashed line (a) = 
absolute number of events) 

b,d: corresponding limits of the mixing angles ^ and $ from the limits of: (b) A eT for several values of Am 2 , (d) A ep " , 
for Am 2 = 2eV 2 /c 4 . For comparison, the exclusion limits from i* e disappearance search at Bugey ^| are shown, also for 
Am 2 = 2eV 2 /c 4 . In b(d) areas to the left(right) of the curves, respectively, and between the horizontal lines are excluded. In 
(d) curves 1 and 2 are the KARMEN limits from — > v & and v^— appearance search. 
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FIG. 13. 90% CL limits to the mixing angles $ and \l/ fixed at Am 2 = 2eV 2 /c 4 deduced from this experiment from the 
oscillation channels Vu,— > v e appearance and i/ e —> v T disappearance in comparison with other experiments: FNAL E531 

v^^Vr appearance |18|, BNL E776 u^— > v e appearance K2| and a 90% likelihood favoured region from LSND Q (shaded 
region) . 
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